Thermal Issues in Casimir Forces Between Conductors and Semiconductors 
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The Casimir effect between metal surfaces has now been well-verified at the few-percent level 
experimentally. However, the temperature dependence has never been observed in the laboratory, 
since all experiments are conducted at room temperature. The temperature dependence for the 
related Casimir-Polder force between an atom and a bulk material has, in contrast, been observed 
between a BEC and a silica substrate, with the environment and the silica held at different tem- 
peratures. There is a controversy about the temperature dependence for the force between metals, 
having to do with the magnitude of the linear temperature term for both low and high temperature, 
the latter being most prominent at large distances. There are also related anomalies pertaining to 
semiconductors. The status of this controversy, and of the relevant experiments, are reviewed in 
this report. 
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I. INTRODUCTION 

The Casimir effect, reflecting quantum vacuum fluctu- 
ations in the electromagnetic field in a region with mate- 
rial boundaries, has been studied both theoretically and 
experimentally since 1948 The forces between dielec- 
tric and metallic surfaces both plane and curved have 
been measured at the 10 to 1 percent level in a vari- 
ety of room-temperature experiments, and remarkable 
agreement with the zero-temperature theory has been 
achieved. For reviews see [2n3l- 

In fitting the data various corrections due to surface 
roughness, patch potentials, curvature, and temperature 
have been incorporated. It is the temperature correction 
that is the subject of the present paper. Temperature de- 
pendence has been detected for the Casimir-Polder force 
between atoms in a Bose- Einstein condensate JfJ @ , the 
theory of which was worked out in Q. A recent experi- 
ment by the Yale group at large distance (0.7-7 /im) Q 
shows the theoretically expected reduction of the high- 
T effect for the Casimir force between metallic surfaces. 
See also @. 

Theoretically, there are subtle issues concerning ther- 
modynamics and electrodynamics which have resulted 
in disparate predictions concerning the nature of these 
corrections. (An overview is given in ficj|.) However, a 
general consensus has emerged that the low-temperature 
correction to the Casimir effect is relatively large, and 
that the linear high-temperature effect should be reduced 
from the naive expectation by a factor of 1/2, in accord 
with the recent Yale experiment (For a critique of 
this experiment, see (TT|-) 
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FIG. 1: Geometry of two semi-infinite parallel dielectric slabs, 
1 and 2, separated by a third dielectric slab of thickness a. 
Each medium is characterized by a permittivity, ej, which is 
a function of the frequency. 



II. CONVENTIONAL APPROACH 

The zero-temperature Casimir effect between parallel 
conducting plates, or between parallel dielectrics, is very 
well understood, and is not controversial. The formula 
for the latter, which includes the former as a singular 
li mit, may be derived by a multitude of formalisms 0- 
HL E Eli- For a system of parallel isotropic dielectric 
media, as shown in Fig. [TJ characterized by a permittivity 
( £ = eH) 




(2.1) 



the Lifshitz force per unit area on one of the surfaces is 
at zero temperature 
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where £ is the imaginary frequency, C = —iw, and the 
longitudinal wavenumber is 
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while the transverse electric (TE) and transverse mag- 
netic (TM) Green's functions are characterized by the 
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respectively, where k\ = Ki/si. 

The attractive Casimir pressure between parallel per- 
fectly conducting planes separated by a vacuum space of 
thickness a is obtained by setting —> oo and £3 = 1. 
In that case the TE and TM contributions arc equal, and 
we have 
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which is Casimir's celebrated result 

The controversy surrounds the question of how to in- 
corporate thermal corrections into the latter result. At 
first glance, the procedure to do this seems straightfor- 
ward. It is well-known that thermal Green's functions 
must be periodic in imaginary time, with period (3 = 1/T 
[Til ]. This implies a Fourier series decomposition, rather 
than a Fourier transform, where in place of the imaginary 
frequency integral above we have a sum over Matsubara 
frequencies 
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the prime being an instruction to count the m = term 
in the sum with half weight. This prescription leads to 
the following formula for the Casimir pressure between 



perfect conductors at temperature T, (t 
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From this it is straightforward to find the high and low 
temperature limits, 
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These are the results found by Lifshitz [15J, Fierz [1 61 ] . 
Sauer [13] , and Mehra [l^l ■ The two limits are connected 
by the duality symmetry found by Brown and Maclay 
[ijj. The pressure may be obtained by differentiating 
the free energy, 
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which takes the following form for low temperature (now 
omitting the exponentially small terms) 
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from which the entropy follows, 



aT<l, (2.10) 



S — 8f F 



3C(3) , 
2tt " 



— T 3 a, aT«l, (2.11) 



which vanishes as T goes to zero, in accordance with the 
third law of thermodynamics, the Nernst heat theorem. 



III. EXCLUSION OF TE ZERO MODE 

However, there is something peculiar about the proce- 
dure adopted above for a perfect metal [20]. It has to 
do with the transverse electric mode of zero frequency, 
which we shall refer to as the TE zero mode. If we ex- 
amine the zero frequency behavior of the reflection coef- 
ficients for a dielectric we see that providing ( 2 e(i() —> 
as C 0, the longitudinal wavenumber m — > k as C 0, 
and hence d — > 00 as £ — > 0. This means that there 
is no TE zero mode for a dielectric. This statement is 
not controversial. However, if a metal is modeled as the 
e — > 00 limit of a dielectric, the same conclusion would 
apply. Because that would spoil the concordance with 
the third law noted in the previous section, the prescrip- 
tion was promulgated fl3l. Il5l| that the e — > 00 limit be 
taken before the ( — ^ limit. But, of course, a real metal 
is not described by such a mathematical limit, so we must 
examine the physics carefully. 

A simple model for the dielectric function is the plasma 
dispersion relation, 
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(3.1) 



where w p is the plasma frequency. For this dispersion re- 
lation, the condition £ 2 e(iC) — > fails to hold as £ — > 0, 
and the idealized prescription result, namely the contri- 
bution of the TE zero mode, follows. 

However, real metals are not well described by this 
dispersion relation. Rather, there is no TE zero mode in 
the Drude model, 
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where the relaxation frequency 7 represents dissipation. 
The Drude model very accurately fits optical experimen- 
tal data for the permittivity for £ < 2 x 10 15 rad/s 
[Hi HI]. For gold, appropriate values of the parameters 
are uj p = 9.03 eV, 7 = 0.0345 eV. 

Let us review the argument by writing the Lifshitz for- 
mula at finite temperature in the form 
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where the second equality uses the Euler-Maclaurin sum 
formula, in terms of 
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where we assume that vacuum separates the two plates 
so K 3 (( m ) = K(Cm) = + Cm- Here the denominators 
are functions of Cm- By changing the integration vari- 
able from m to Cm we immediately see that the integral 
term in the Euler-Maclaurin sum formula corresponds 
precisely to the zero-temperature result. 

One must, however, be careful in computing the low 
temperature corrections to this. One cannot directly ex- 
pand the denominator d in powers of C because the k± 
integral ranges down to zero. Let us rewrite the TE term 
there as follows: 
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Evidently, for the Drude model, or more generally, when- 
ever 



lim C 2 [e«) - 1] = 0, 
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the TE zero mode vanishes, /( TE )(0) = 0. However, it 
is important to appreciate the physical discontinuity be- 
tween m = and m = 1 for room temperature. At 300 
K, while Co = 0, Ci = 27rT = 0.16 cV, large compared 
the relaxation frequency 7. Therefore, for m > 0, 
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provided the significant values of Cm and k are small 
compared to the plasma frequency lu p . This is just the 
ideal metal result. Insofar as this is accurate, this ex- 
pression yields the low- and high-temperature corrections 
seen above. However, there is now a discontinuity in the 
function /( TE ). As Cm -> 0, 
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rather than zero. 

This implies an additional linear term in the pressure 
at low temperatures: 
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Exclusion of the TE zero mode will also reduce the linear 
temperature dependence expected at high temperatures, 
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one-half the usual ideal metal result seen in Eq. (I2.8al) . 



IV. PROXIMITY FORCE APPROXIMATION 

Most experiments are carried out between a sphere (of 
radius R) and a plane. In this circumstance, if R a, 
a being the separation between the sphere and the plate 
at the closest point, the force may be obtained from the 
proximity force approximation (PEA) (23|, 



& = 2irRF(a), 



(4.1) 



F(a) being the free energy for the case of parallel plates 
separated by a distance a. Thus in the idealized descrip- 
tion, the low temperature dependence including our lin- 
ear term is for aT <C 1, 
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Since this conversion is trivial, in the following we will 
restrict attention to the straightforward parallel plate sit- 
uation. (Recently, use of PFA has been superseded by nu- 
merical/analytical methods, e.g., multiple-scattering and 
boundary element techniques. See [2~H for a review.) 

These results are only approximate, because they as- 
sume the metal is ideal except for the exclusion of the 
TE zero. Elsewhere, we have referred to this model as 
the Modified Ideal Metal (MIM) model Evidently, 
for sufficiently low temperatures the approximation used 
here, that Ci 3* 7> breaks down, the function f(m) be- 
comes continuous, and the linear term disappears. In- 
deed, numerical calculations based on real optical data 
for the permittivity show this transition. Because this 
linear behavior does not persist at arbitrarily small tem- 
peratures, it is clear that the conflict with the third law 
anticipated in the arguments in the previous section do 
not apply. In fact, as we shall see in Sec. IVI| the en- 
tropy does go to zero at zero temperature for any realis- 
tic metal. Only in special idealized model cases does the 
entropy problem persist, wherein Ci S> 7 remains true all 
the way to zero temperature; i.e., when 7 tends to zero 
as 7 cx T or faster as T -> [H, S3- 



V. KRAMERS-KRONIG RELATION 

As noted above, there are strong thermodynamic and 
electrodynamic arguments in favor of the exclusion of 
the TE zero mode. Essentially, the point is that a re- 
alistic physical system can have only one state of low- 
est energy. Electrodynamically, one can start from the 
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Kramers-Kronig relation that relates the real and imag- 
inary parts of the permittivity, required by causality, 
which can be written in the form of a dispersion rela- 
tion for the electric susceptibility [28j 



xM 



v 
-in 



du) - 



(5.1) 



If the spectral function p(u>') > is nonsingular at the 
origin, it is easily seen that w 2 x(u;) — > as cj — > 0, which 
as shown in Sec. IIIII implies the absence of the TE zero 
mode. Conversely, p(u>') must have a <5-function singu- 
larity at the origin to negate this conclusion. This would 
seem implausible for any but an overly idealized model. 
In contrast, in the Drude model the spectral function 
is non-singular, becoming singular only in the plasma- 
model limit: 
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As noted above, use of the plasma model in the reflec- 
tion coefficients would lead to the conventional temper- 
ature dependence, but this dispersion relation is incon- 
sistent with real data. It may be argued (2(| that in the 
ideal Bloch-Griineisen model [30[ the relaxation param- 
eter goes to zero as 7 oc T 5 near zero temperature, thus 
causing the entropy to tend to a non-zero value at zero 
temperature. Physically this idealized case was shown 
to correspond to a glass-like state of frozen, randomly 
oriented bulk currents, whose ground state is strongly 
degenerate, hence retaining finite entropy in this limit 
without thereby violating the third law [3l|, HH . 

More to the point, however, real materials exhibit scat- 
tering by impurities and certainly from boundaries, 1 so 
that 7 remains finite. In any case, at sufficiently low tem- 
peratures the residual value of the relaxation parameter 
does not play a role, as the frequency characteristic of 
the anomalous skin effect becomes dominant [34[ . More- 
over, the Purdue group [35j also extrapolate the plasma 
formula from the infrared region down to zero frequency, 
whereas in fact frequencies very small compared to the 
frequency corresponding to the separation distance play 
a dominant role in the temperature dependence. 

Finally, we emphasize that all present Casimir exper- 
iments are carried out at room temperature, where the 
known room temperature data are relevant. 



VI. THIRD LAW OF THERMODYNAMICS 

The principal reason for the theoretical controversy has 
to do with the purported violation of the third law of 



1 Indeed, for very pure metals 7 is found experimentally to be 
sample-size dependent, and hence reaches a nonzero value even 
for an imperfection-frce metal when its Bloch-Griineisen electron 
mean free path becomes comparable to sample dimensions 33]- 
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FIG. 2: Numerical evaluation of the free energy between two 
gold halfspaces as a function of temperature. The inset gives 
details for low T. 



thermodynamics if the TE zero mode is not included. 
If ideal metal reflection coefficients are used otherwise 
(the MIM model) such a violation indeed occurs, because 
the free energy per unit area for small temperature then 
behaves like 



F = F +T 
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However, we and others have shown that for real metals, 
the free energy per area has a vanishing slope at the 
origin. Indeed, in the Drude model we have Q 
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for sufficiently low temperatures. The T 5 / 2 correction 
to this has been calculated in [36[ . Figure [2j from that 
reference, shows a numerical evaluation of the free energy 
for gold, and how the low-temperature linear behavior 
gives way to a flat slope for very small temperatures. The 
exact mathematical requirements for the formal violation 
of the third law are reviewed in [37| . 

There is, however, an intermediate range of tempera- 
tures where it is expected that the entropy is negative. 
We do not believe that this presents a thermodynamic 
difficulty, and reflects the fact that the electrodynamic 
fluctuations being considered represent only part of the 
complete physical system [HI . It was also pointed out by 
Pitaevskii that the "Casimir entropy" is not a measure 
of the system's absolute entropy, but is to be understood 
as the difference between the (positive) entropies of an 
initial and final state, and there is no theorem dictating 
the sign of such a difference [3!| . 

Independent arguments also lend support to our point 
of view. Jancovici and Samaj [4p| ] and Buenzli and 
Martin [4l], [42[ have examined the Casimir force be- 
tween ideal-conductor walls with emphasis on the high- 
temperature limit. Not surprisingly, ideal inert boundary 
conditions arc shown to be inadequate, and fluctuations 
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within the walls, modeled by the classical Debye-Hiickcl 
theory, determine the high temperature behavior. The 
linear in temperature behavior of the Casimir force is 
found to be reduced by a factor of two from the behav- 
ior predicted by an ideal metal, just as seen above. This 
is precisely the signal of the omission of the m = TE 
mode. Thus, it is very hard to see how the correspond- 
ing modification of the low-temperature behavior can be 
avoided. 

It is just the reduction of this dependence by the factor 
of 1/2 that is seen in the Yale experiment Q. 



VII. BOHR- VAN LEEUWEN THEOREM 

Bimonte [43| has suggested that the Bohr- van Leeuwen 
theorem, which says that classically at thermal equilib- 
rium electromagnetic fields decouple from matter, ex- 
plaining why normal metals do not exhibit strong dia- 
magnctism, requires that the TE reflection coefficient 
vanishes as zero frequency. Plasma-like prescriptions are 
thus ruled out. The Bohr-van Leeuwen theorem does not 
apply to magnetic materials or superconductors, since 
quantum effects dominate there. Thus experiments on 
these materials might reveal quite different behaviors. 
Nevertheless, the Bohr-van Leeuwen theorem would seem 
a more natural criterion to turn to than the third law, 
for systems where magnetic effects are negligible at room 
temperature; the former applies in the high temperature 
limit which is exactly where the two Casimir predictions 
differ the most, whereas the latter applies only as the 
temperature tends to zero, a situation very far from lab- 
oratory conditions in any experiment reported to date. 



Casimir force between two Al half-spaces 
relative to ideal result 




FIG. 3: Temperature dependence of the Casimir force be- 
tween aluminum plates. The abscissa is the separation in 
nanometers. 
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VIII. EXPERIMENTAL CONSTRAINTS 

We have marshalled theoretical arguments that seem 
to us quite overwhelming in favor of the absence of the TE 
zero mode in the temperature dependence of the Casimir 
force between real metal plates, which seem to imply un- 
ambiguously that there should be large (~ 15%) thermal 
corrections to the Casimir force at separations of order 1 
micrometer. Calculations based on this theory, and using 
real optical data for aluminum, for example, have been 
given, as shown in Fig.[3]for aluminum plates, taken from 
[10j . The difficulty is that, experimentally, it is not easy 
to perform Casimir force measurements at other than 
room temperature, so current constraints on the theory 
all come from room temperature experiments. Then all 
one can do is compare the theory at room temperature 
with the experimental results, which must be corrected 
for a variety of effects, such as surface roughness, finite 
conductivity, and patch potentials. A deviation between 
the zero temperature theory and room temperature ob- 
servations is a measure of the temperature correction. 

Figure 3] shows the expected behavior of the pressure 
as a function of plate separation a at room temperature. 



FIG. 4: The pressure between two gold plates multiplied a 4 , 
plotted versus the plate separation a for room temperature. 
Data for the permittivity is taken from [2l|, [2^] . This figure 
is taken from [25| . 



Since, approximately, the dependence is through the di- 
mensionless factor (in natural units) aT, the high temper- 
ature linear behavior is translated into a linear behavior 
in the separation. (The short distance linear regime is 
also seen with the expected slope.) 

The temperature correction is evidently relatively 
largest at the largest separations, where, unfortunately, 
the total Casimir force is weakest. Lamoreaux's new ex- 
periment [|[ is conducted at the 0.7-7 /mi scale, which 
is where the linear high-temperature result is dominant. 
At these large distances, the Casimir force is very small, 
and is overwhelmed by patch potentials. Nevertheless, af- 
ter the patch-potential effect is extracted and subtracted 
from the data, the effect is to strongly favor the Drudc- 
model, that is, exclusion of the TE zero mode, as illus- 
trated in Fig. [4] The procedures and conclusions of Q 
are strongly critized in [Til ] . 

It is the most precise experiments of the Purdue group 
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[4:4 . |45| that claim the extraordinarily high precision to 
be able to see our effect at distances as small as 160 
nm. Indeed, they see no deviation from the corrected 
zero-temperature Lifshitz theory using optical permittiv- 
ities, and hence they assert that our theory (the "Drudc 
model") is decisively ruled out. The effect we predict for 
the temperature correction is only at the 1% level at a dis- 
tance of 160 nm, so the measurement must be performed 
to that accuracy to sec the effect there. (For the usually 
employed sphere-plate configuration, AJ^ « 2.5% at 
a = 160 nm.) Taking the reported accuracy at face value, 
the nondissipative plasma-model approach is favored at 
a very high confidence level. While doubt has been aired 
on a general basis that a 1% accuracy can be possible 
phi , we are not aware of any criticism regarding concrete 
points in Decca's impressive experiment, and the results 
found must be regarded as something of a mystery. 

Recently, there has also been a report from the River- 
side group [13, claiming that the Drude model is also 
disfavored, although not at as high a level of significance. 
We cannot comment further on this result at present, ex- 
cept to note that the theoretical support for both the 
Purdue and the Riverside experiments is identical. 



IX. ANOMALY FOR SEMICONDUCTORS 

Mostcpancnko and Klimchitskaya have also claimed 
there is a similar anomaly, now affecting the TM reflec- 
tion coefficient, for dielectrics and semiconductors 48|. 
Here is a simple way to understand that argument 49J. 
Suppose we model a dielectric with some small conduc- 
tivity by the permittivity function 
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(The only essential point is that as £ — > 0, e — > eq if a = 
0, otherwise e — > oo.) The Casimir (Lifshitz) pressure is 
given by 
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For the case of an ideal metal, it was r TE which was 
discontinuous: 



r TE (C = 0)=0, 
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so this gave a linear temperature term when the sum over 
Matsubara frequencies is converted to an integral accord- 
ing to the Euler-Maclaurin formula, for example. For a 
dielectric the TE reflection coefficient is continuous, but 
if there is vanishingly small (but not zero) conductivity 
the TM coefficient is discontinuous: 
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This leads to a linear T term in the pressure: 
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where the logarithmic integral function is 
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Note that the linear term vanishes for £o —> oo. The free 
energy term F is obtained from this by multiplying by 
a/2. Thus at zero temperature, the entropy is nonzero, 
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as found in [48| , which again appears to violate the third 
law of thermodynamics, stating that the entropy of a 
nondegencrate system must vanish at zero temperature. 

We believe that real materials cannot undergo this vio- 
lation, and that, properly treated, the entropy must van- 
ish at zero temperature. The details of this assertion, 
however, have not been worked out to our knowledge. In 
materials with a very small concentration of free charges, 
it seems highly unlikely that a simple, local, Drudc-typc 
model can provide an adequate description of reflectivity 
[2r| , and salvation may well be found to lie in a complete 
inclusion of spatial dispersion effects. 

Three simultaneous and very similar efforts were made 
in this direction j50Tj52| . These partly seek to bridge the 
gap between the predictions of a naive Drude-likc model 
(19.111 an d an experiment performed by Mohideen's group 
[53l |54| . In that experiment, the carrier concentration 
of a semiconductor was increased by laser illumination, 
and the observed behaviour of the Casimir force is sur- 
prising at first sight: 1) When the carrier concentration 
is high, the conductivity must be included if the sim- 
ple model (|9.1|) is used. 2) When the carrier density is 
low (but nonzero), the conductivity must be excluded in 
order to reproduce measured data. In light of the spa- 
tial dispersion theory based on Debyc-Hiickel screening 
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[50U521 |55| . such behavior is not so surprising, however, 
although these theories fail to reproduce the data per- 
fectly. Unlike the mystery of Decca's experiment with 
metal surfaces, the experimental result here appears intu- 
itively reasonable, but a complete theoretical description 
without any ad hoc prescriptions is still lacking. 

X. REAL TEMPERATURE EXPERIMENTS 

In any case, it would seem imperative to perform 
experiments at different temperatures in order to pro- 
vide evidence for or against temperature dependence of 
Casimir forces. However, such experiments are extremely 
difficult, and, for example, Onofrio has abandoned his 
attempt to do such a measurement (private communica- 
tion). For example, the increase of temperature from 300 
K to 350 K should reduce the force by about 2.5%, for 
Au-Au plates when the separation is 1.0 micrometer. 

Even though the new Yale experiment Q is consistent 
with the revised theory, it needs to be independently con- 



firmed, especially in view of the very large electrostatic 
corrections that must be supplied. Also, direct measure- 
ment of the temperature dependence would be highly 
beneficial. (The only direct temperature measurement 
of temperature dependence is that for the CP force on 
a BEC of Rb atoms 0.) We encourage experimentalists 
to redouble their efforts, for the issues involved touch at 
the heart of our fundamental theoretical understanding 
of electrodynamics, statistical mechanics, and quantum 
field theory. 
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